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We evaluate quantum corrections to conductivity in an electrically gated thin film of a three- 
dimensional (3D) topological insulator (TI). We derive approximate analytical expressions for the 
low-field magnetoresistance as a function of bulk doping and bulk-surface tunneling rate. Our results 
reveal parameter regimes for both weak localization and weak antilocalization, and include diffusive 
Weyl semimetals as a special case. 
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I. INTRODUCTION AND OVERVIEW 

The theoretical discovery^ of 3D topological insulators 
(TIs) in 2006 precipitated an avalanche of experiments 
aimed at detecting the signature behavior of these uncon- 
ventional solids. Since then, angle-resolved photoemis- 
sion spectral have given evidence for the Dirac-like dis- 
persion and the momentum-dependent spin texture of TI 
surface states, whereas local STM probes have indicated 
a characteristic suppression of backscattering off surface 
imperfections^ However, the most desired observation of 
a hallmark dc conduction confined to the surface layer of 
a 3D TI remains elusive.— The main problem is conduc- 
tion through the bulk: 3D TIs are narrow-gap semicon- 
ductors, rich in bulk carriers that are either thermally 
activated and/or donated by crystalline lattice imperfec- 
tions. Along with attempts to reduce bulk charge carri- 
ers, experimentalists are developing techniques which al- 
low to register a separate conduction channel on the sur- 
face of a 3D Tl£ Chief among these are measurements of 
low-field magnetoresistance combined with electrostatic 
gating of thin-film samples 

Low-field magnetoresistance measurements unveil the 
interference correction 8a to the Drudc conductivity 
od<^ At low temperatures, ao is defined by indepen- 
dent acts of scattering of electrons off the crystal's im- 
perfections, and is proportional to the classical diffu- 
sion constant D. When the phase relaxation length 1^ 
is parametrically longer than the scattering mean free 
path, quantum interference affects the conductivity to 
a measurable extent. The sign of the interference cor- 
rection depends on the strength of spin-orbit interac- 
tions. For weak spin-orbit interactions (l so 3> 1$, where 
l so is the spin-orbit scattering length), it follows that 
8a < 0. This is called weak localization (WL). In con- 
trast, strong spin-orbit interaction (l so <§C 1$) leads to 
suppression of backscattering and thus 8a > 0. This is 
called weak antilocalization (WAL). Being interference 
effects, WL and WAL are degraded by a magnetic field 
H when H > = $ /(8nlV), where $ = h/e is the 
flux quantum. Yet, ctd is nearly immune to H at such 
low fields. Therefore, the low-field magnctoconductivity 
reads Aa(H) = a{H) - a(0) ~ 8a(H) - 8a(0). 



All experiments to date report WAL in 3D TI thin 
films' and ascribe it to the strong spin-orbit interac- 
tion in the electronic bands of these materials. For film 
thickness less than 1^, the measured Aa(H) agrees well 
with the functional form provided by 2D WAL theory, 
namely 



Aa{H)~a{e 2 /2K 2 h)f{H !l> /H), 



(1.1) 



where f{z) = lnz — '0(1/2 + z), with tp and a being the 
digamma function and a number^ respectively. In a 
system with a single conduction channel, a is universal 
and equals 1/2. The WAL contributions add for systems 
which are isolated from each other. For example, having 
two independent parallel conduction channels yields a = 
1, irrespective of the ratio of Drude conductivities of the 
two subsystems. 

The relation between a and the number of paral- 
lel channels is at the heart of recent magnetoresistance 
experiments in 3D TIs ^ 10 ! 11 Overall, the coefficient a 
is found to depend on the gate voltage. For some 
devices ) 7 ' 10 : 11 it changes from a = 1/2 all the way to 
a = 1. A plausible interpretation for this variation is 
presented in Ref . [ll[ . At zero or positive bias applied to 
the top gate, electrons from the n-doped bulk reach the 
surface states easily: the entire film acts as a single elec- 
tron system, and a = 1/2. At negative bias, electrons are 
repeled from the top surface and, for strong enough bias, 
a depletion layer is formed adjacent to it. This deple- 
tion region separates the film into two subsystems: bulk 
carriers (combined with surface carriers from the bottom 
surface) on one side, and top-surface carriers on the other 
side. For a wide enough depletion layer, a = 1. 

In spite of the ongoing scrutiny on the experimental 
front, quantum corrections to conductivity in 3D TIs 
have stimulated relatively little theoretical activity. Even 
though the WAL contribution from TI surface states 
has been calculated explicitly ; 16 ! 17 there are no calcula- 
tions that incorporate conducting 3D bulk states. The 
main reason for this omission may be the prevailing 
view that quantum corrections originating from bulk TI 
states ought to be conceptually identical to those in or- 
dinary strongly spin-orbit coupled systems, i.e. of WAL 
type. Recently, an objection to this viewpoint has been 
raised^ declaring that quantum well states in ultrathin 
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TI films may contribute via WL rather than WAL. Al- 
though suggestive, the calculation of Ref. [lU is limited 
to quasi- 2D films and disregards the coupling between 
bulk and surface states, which leaves out several experi- 
ments of interest. Besides, its extrapolation to 3D bulk 
states has not been carried out properly. 

In this paper we evaluate Act for gated 3D thin films, as 
a function of the bulk carrier concentration and account- 
ing for the coupling between surface and bulk states. Our 
calculation applies to TI films that are thicker than the 
bulk mean free path, thinner than 1$, and not highly 
doped. In these films, bulk carriers are three-dimensional 
and are concentrated around the T point of the electronic 
band structure. The resulting approximate analytical 
expressions for Act (Eqs. (|235|) . ([XT?]) and ([3~T5]) ) are 
aimed at improving the interpretation of magnetoresis- 
tance measurements in TIs, in Weyl semimetals^ and 
in some class of topologically trivial materials. Although 
a few of our observations resemble those developed for 
graphene^ and 2D TIs,— there are qualitative differences 
originating from the 3D Dirac nature of bulk carriers in 
3D TIs. 

Altogether, the results reported here paint a richer pic- 
ture than previously anticipated. On one hand, we con- 
firm the conventional crossover between a = 1/2 and 
a = 1 as a function of the gate voltage: the former corre- 
sponds to the case of coherently-coupled bulk and surface 
electron states, while the latter indicates a single decou- 
pled Dirac cone on the top surface along with generic 
WAL from the rest of the film (containing coupled bulk 
and bottom surface). On the other hand, less conven- 
tional results arise when the Fermi energy is close to the 
bulk band edge or when the Fermi energy is much larger 
than the bulk bandgap: in the former regime the bulk 
exhibits WL with a = — 1, whereas in the latter regime 
the bulk exhibits an anomalous WAL with a = 1. These 
two "unusual" bulk regimes, combined with the surface 
contributions, may result in a range of a including a < 
and a > 1. 

The rest of this work is organized as follows. In Sec- 
tion II we evaluate quantum corrections to bulk conduc- 
tivity. Readers not interested in technical details should 
read subsection IIA and quickly scan through IIB and 
IIC in order to get acquainted with the nomenclature; 
the main results of the section are collected in Section 
IID. The well-known message from IIA is that at low en- 
ergies bulk electrons of TI films behave as massive 3D 
Dirac fermions with spin and valley (or orbital) degrees 
of freedom. The direction of spin is locked with that of 
momentum, and valleys are coupled to one another by the 
mass of the Dirac fermions. The special case in which the 
Dirac mass vanishes is a time- and inversion-symmetric 
Weyl semimetal. 

In Section IID we identify and count the number of 
"soft" Cooperon modes, which determine the magnitude 
and sign of Act in the bulk. Each soft Cooperon obeys 
a classical difussion equation and is thus associated with 
a conserved physical quantity. Since charge is conserved, 



there is at least one soft Cooperon in (non-magnetic) bulk 
TIs. We find that additional soft Cooperons can emerge 
depending on the bulk doping concentration as well as 
the bulk bandgap. This realization leads to the most im- 
portant results in IID, Eqs. (|2.33[) - (|2.35[) . which indicate 
that for bulk states a may acquire three different univer- 
sal values. On one hand, WL with a = — 1 is possible 
when the bulk Fermi surface is "small" (as defined in the 
text), because in this case the spin- momentum locking 
of bulk states becomes weak and the spin of electrons 
is nearly conserved. In contrast, WAL with a = 1 can 
arise for bulk TIs with particularly small bandgaps, be- 
cause in such case bulk electrons can be described by a 
3D analogue of graphene with two nearly decoupled val- 
leys, each contributing 1/2 to a. For a more generic case, 
in which neither valley nor spin arc approximately con- 
served, the quantum interference is similar to that of an 
ordinary film with strong spin-orbit coupling and there- 
fore a = 1/2. Magnetic fields perpendicular to the TI 
film can be used to induce crossovers between different 
universal regimes of a. The accessible values of a and the 
corresponding crossover fields depend on the bulk elec- 
tron density. 

In Section III we evaluate the full Act in 3D TI thin 
films, which comprises coupled bulk and surface contri- 
butions. Sections III A and IIIB cover preliminary ma- 
terial that is needed to derive the main results in IIIC. 
Section III A reviews the well-established fact that, in ab- 
sence of magnetic order, isolated TI surface states exhibit 
WAL with a = 1/2 (in this paper we assume one Dirac 
cone per surface). Section IIIB develops a diagrammatic 
framework for evaluating quantum corrections to conduc- 
tivity in ordinary tunncl-couplcd layers. Readers who 
arc not interested in technicalities can disregard the dia- 
grams in the figures and concentrate on the outcome of 
the calculation (Eqs. (|3.8[) - (|3.12[) 'l. as well as on the sub- 
sequent discussion. One qualitative point made therein 
is that the crossover from weak to strong coupling (which 
is accompanied by a change in a from 1 to 1/2) occurs 
when the interlayer resistance for a square of area l^ be- 
comes smaller than the sum of the classical intralaycr 
resistances. 

Section IIIC combines results from IID, IIIA and IIIB 
in order to figure out quantum corrections to conductiv- 
ity in experimentally realized TI films. The most impor- 
tant results in IIIC are Eqs. (|3.17[) and (|3.19[) . which 
describe how Act depends on the bulk doping concen- 
tration, on the phase relaxation rate, and on the bulk- 
surface tunneling rate. Some special cases of these re- 
sults arc highlighted in Appendix [F] A salient conclusion 
is that the WL regime of isolated bulk states is gener- 
ally eliminated when either one of the film surfaces is 
strongly coupled to bulk states, in which case the film 
displays 1/2 < a < 1. However, WL can still be present 
if the TI surfaces have short phase relaxation lengths. 

Finally, Section HID characterizes the electrostatics of 
the depletion layer and estimates the bulk-surface tun- 
neling rate in TI films. This estimate confirms cxperi- 
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mental indications showing that both weak and strong 
bulk-surface coupling are accessible by mediation of a 
gate voltage. 



II. QUANTUM CORRECTIONS TO BULK 
CONDUCTIVITY 

This section is devoted to evaluating 5a for the bulk 
states of a 3D TI. As a byproduct, wc derive 5a for a time- 
reversal symmetric Weyl semimetal. The contribution 
from TI surface states will be discarded until the next 
section. 




FIG. 1: Bulk energy bands of an n-doped 3D TI near the 
r point, in the spherical approximation. The momentum k 
is measured from the V point. The energies 6f and M are 
measured with respect to midgap. 



A. Model 

The bulk band structure of a 3D TI near the T point 
can be approximated by the following kp Hamiltonian^i 

k 

/i(k) ~ e(k)l 4 + M(k)l 2 t z + H (v z k z a z + v x k x • tr 1 -) t x , 

(2.1) 

where r is an orbital pseudospin (t z = T, £>), <x is the real 
spin (a z =j, J,), k = (kj_, k z ) is the momentum measured 
from the T point of the Brillouin zone, In is an N x 
N identity matrix, * = (* Tt , $ TJ _, \& Bt , ^ Bi ) is a 4- 
spinor, e(k) = e(— k) is the part of the Hamiltonian that 
is independent of spin/pscudospin indices, v z and v± arc 
the Fermi velocities, and M(k) = M — M%k\ — M 2 k 2 is 
the mass term (independent of spin). Mo, M\ and M 2 
are constants. 

Equation (|2.1j) captures the bottom of the conduction 
band and the top of the valence band in the vicinity of 
the r point (k = 0), where the bandgap is smallest. It 
models 3D Dirac fcrmions with a Dirac mass that equals 
half the energy gap. For the purposes of this paper we 
ignore e(k), and assume M(k) = M = const > as well 
as spherical symmetry (v z = vj_ = v). These assump- 
tions simplify calculations without incurring in qualita- 
tive loss of generality. For instance, the XXZ anisotropy 
can be modeled by promoting the diffusion constant from 
a scalar to a matrix. Also, the k 2 terms in M(k) can 
be incorporated into our final results by M -+ |M(kj?)|, 
where kp is the Fermi wave vector. Note that in absence 
of spherical symmetry the Fermi surface does not have a 
constant mass; this complication will be disregarded in 
the present paper. Finally, e(k) can be absorbed into the 
definition of the Fermi energy. 

The energy eigenvalues for ft.(k) in the spherical ap- 
proximation arc Ek± = ±\Jh 2 v 2 k 2 + M 2 , each doubly 
degenerate (Fig. [J). The corresponding Bloch states can 
be written as 



where V is the volume of the TI and a <G {1,2,3,4} is 
a band index (1 and 2 denote conduction bands, while 3 
and 4 denote valence bands). This a is obviously unre- 
lated to that of Eq. (jl.ljl : from here on it will be clear 
from the context which one we are referring to. For con- 
creteness we set the chemical potential in the bulk con- 
duction band, although all results obtained below will 
be directly applicable to p-doped bulk TIs as well. The 
density of conduction band electrons is then 



M 2 ) 



3/2 



TT 2 h 3 V 3 



(2.3) 



where cf is the Fermi energy measured from the middle 
of the bulk energy gap. Adopting the basis {\T f), |T J, 
),\B f), \B I)}, the two eigenspinors corresponding to 
the conduction bands near the T point are 



Uk) = ,/^±^fi,o, hvk+ 



2E k 



E k + M ' E k +M 



12k) 



Ek + M ( hvk- -hvk z 

0,1, , - 

2E k \ Ek + M 7 E k + M 



,(2-4) 



(1/VV) exp(ik-r)|ak), 



(2.2) 



where k± — k x ± ik y and Ek = ^k,+- Since all non-Hall 
dc transport properties of good conductors involve states 
close to the Fermi energy, we hereafter ignore valence 
bands. 

Unlike in the k • p Hamiltonians for graphene and 2D 
(or quasi- 2D) TIs, Eq. (|2.ip cannot be decomposed into 
two 2x2 block diagonal matrices (due to M ^ 0). In 
addition, the k z band dispersion absent in 2D cannot be 
neglected in our case. These two features make the cal- 
culations and results of this section quite different from 
those of Refs. [l^Iilll^. 

Equation (|2.1j) becomes inaccurate when the chemical 
potential moves up in the conduction band and electron 
pockets away from the T point begin to be populated. 
These additional pockets contribute to quantum inter- 
ference, and the total 5a depends on the scattering rate 
between different electron pockets. Although a realistic 
study of the full band structure is beyond the scope of 
this paper, we expect calculations based on Eq. (|2.1I) to 
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provide a generic understanding of quantum corrections 
to conductivity in 3D Dirac materials at low-to-moderate 
doping concentrations. 



B. Formalism 

In order to quantify the conductivity of a bulk TI, 
we begin by characterizing the simplest possible disorder 
potential: Vdis(r) = V(r)l4, which is time-independent 
(clastic) and independent of spin as well as orbital de- 
grees of freedom. For simplicity we assume V(r) to be 
slowly-varying at the atomic scale, yet short-ranged com- 
pared to the mean free path: V(r) — Vb<5(r). It is due to 
its slow spatial variation on atomic lenghtscales that Vd; s 
becomes an identity operator in orbital space. With such 
disorder realization, the Fermi-surface lifetime To for the 
a = 1, 2 eigenstates in Eq. (|2.2j) obeys 



1 

ro 



2ttuq 

h 

TTUqV 



j ^2\( ak F\ a ' k, F)\ 2 S^F - E Wa ,) 



M 2 



(2.5) 



where L = j d 3 k/(2ir) 3 , uq = n^Vy 2 , n, is the density of 
impurities, and v is the density of states per band and 
per unit volume at e F . 

A related quantity is the transport scattering rate r _1 , 



1 _ 27T"Uo 

r 



j 5^(1 - k F • k^)|(ak F |c/k F )| 2 <S(e F - E k , a ,) 



2 e 2 F + 2M 2 
3^ e 2 F + M 2 ' 



(2.6) 



The momentum-dependence of |ak) makes to ^ r even 
for (5-function impurity potentials. Throughout this work 
we impose (e F — M)t 3> h or equivalently kpl 3> 1, where 
/ = (Dt) 1 / 2 is the elastic mean free path, 



k F = (e 2 F - M 2 )^ 2 /(hv) 



is the Fermi wave vector and 



D 



pt/3 = v 2 t{1 - M 2 /e 2 F )/3 



(2.7) 



(2. 



is the classical diffusion constant. 

Next, we consider a TI with spatial dimensions L x L 
in the xy plane and a thickness W along the z direction. 
We take a thin film geometry with L ^> 1$ ^> I and 
(^>F>1, where 1$ = {Dt^) 1 / 2 is the coherence length 
and is the phase relaxation time. The conductivity of 
this film is 



CD 



8a, 



(2.9) 



where ap, is the classical (Drude) part and 8a is the part 
coming from quantum interference. 
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FIG. 2: (a) Feynman diagram for Sai, defined in the text. 
Filled squares denote velocity operators (including disor- 
der vertex corrections), C is the Cooperon. (b) Diagram- 
matic representation of the Bethe-Salpeter equation for the 
Cooperon. Crosses correspond to impurity scattering cen- 
ters. Solid lines with arrows are disorder-averaged Green's 
functions, (c) Additional Feynman diagrams that contribute 
to conductivity of 3D TIs even when impurity scattering is 
isotropic. 



On one hand, the Drude conductivity can be approxi- 
mated as 

a D ^pY, [ ^(k)^ a (k)Gf (k)G#(k), (2.10) 

where we have assumed a spatially uniform dc electric 
field and a,/3£ {1,2}. v*g = (ak|v • x\(3k) is a matrix 
element for the x-component of the bare velocity operator 
v = vt x ct, which obeys 

v Q(9 (k) =hv 2 (k/E k )8 af) (for a,/? e {1,2}). (2.11) 

Disorder vertex corrections renormalize Eq. (|2.11l) to 



v Q/ 3 = v q/j (t/t ), 
see Appendix |A"1 In addition, 



G^ A) (k) 



ih 
2ro" 



~i -l 



(2.12) 



(2.13) 



is the zero-frequency retarded (advanced) Green's func- 
tion in the band eigenstate basis. Using G± (k) = 
Gf A \k) = G R ( A \k), Eq. (I2~10l) yields 



a D = 2e 2 vD. 



(2.14) 



On the other hand, the quantum correction 8a can 
be written as 8a ~ 8a \ + 8ai , represented pictorially in 
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Fig. [2j Following standard approximations, the expres- 
sion for <5<7i is 

Sa ^^ £ /^(k)e|, a ,(-k)Gf(k)G«,(-k) 



and 



1 



=^(k)G^(-k)^/ 7 ^^(k ) k, Q ). 



WJ (2tt) 2 



(2.15) 



In the second line of Eq. (|2.15[) we have exploited the con- 
dition W <C 1$, which allows to set Q z = everywhere. 

Cf^(ki, k 2 , Q) are the matrix elements of the Cooperon 

matrix C in the band eigenstate basis. Q = (Q x ,Qy) is 
the momentum of the Cooperon, whose magnitude ranges 
from to ~ {Dt)- 1 / 2 . C obbeys the Bethe-Salpeter 
equation (Fig. [2p): 



C^(k 1 ,k 2 ,Q)=r^(k 1 ,k 2 ,Q) 



0/3' 



rfC(ki,k3,Q) 



x G£„(k 3 )G£„(-k 3 + Q)Cf:f (k 3 ,k 2 , Q) 



(2.16) 



where a sum over repeated indices is implied and 
r^(k 1; k 2 , Q) = Uo^kxl/3' - k 2 + Q)(o/ - k x + Q|ak 2 ) 

is the bare disorder vertex (first term on the right hand 
side of Fig. Eb). 

Equation (|2.16|) is a complicated integral equation be- 
cause C^f is a function of three momenta. This is un- 
like in simplest examples, where the Cooperon depends 
only on Q. The difficulty originates from the momentum- 
dependence of |ak), which cannot be overlooked as it cru- 
cially determines both the magnitude and the sign of 8a. 
One procedure^ to solve Eq. (|2.16p starts by writing the 
Cooperon in the two-particle spin/orbit basis {\m,m')}, 
where mG {T f, T I, B f , B i} : 

Cjf a Qi 1 M,Q)= J2 («',-ki + Q|m')<)8ki|m} 

mjm', n,n' 

x (n|/3',-k 2 + Q)(n'|ak 2 )C;X'(Q)- (2-17) 

We then make the ansatz that C™Jl, depends on Q but 
not on ki and k.2', the entire ki- and k 2 -dependence of 
C^(ki,k2,Q) originates from the overlap matrix ele- 
ments of Eq. (|2.17[) . The internal consistency of this 
ansatz can be demonstrated by substituting Eq. (|2.17[) 
in Eq. (|2 . 10|) . which produces an algebraic equation for 
C™,™ , that is more tractable than the original integral 
equation: 

Cv(Q) = «o^ + E^i'(QA(Q). ( 2 - 18 ) 



where 



/ii 3 k 



(2.19) 



G^r(k) = ]T(m|ak}G5 A) (k)(ak|0. 



(2.20) 



In matrix language, Eq. (|2.18j) can be rewritten as 

C=(1 W -U)- 1 C^\ (2.21) 



where = Uolxe- 

Eq. (|2.1T[) in order to recover C^, a . During this opera- 
tion we neglect Q in the overlap matrix elements, which 
is a good approximation because 8a is dominated by el- 
ements of C™,™ ,(Q) that are strongly peaked at Q ~ 0. 

For ef in the conduction band, we once again limit 
ourselves to a,/3,a',0' G {1,2} in Eq. (|2.15j) . Then we 
can approximate k ~ k^ inside the Cooperon, and an 
integration over |k| yields 



Sa\ 

where 
C(Q)e 



d 2 Q -r, 



e 2 1 
h 2 W J {2ir) 2 



C(Q), 



(2.22) 



dQw £ 2 

An 



q.q' = 1,2 



Cr a (kF,k^,Q) (2.23) 



and cZfik is the differential solid angle subtended by k. 

Note that 8a\ depends on the lifetime To of Bloch states 
as well as on the transport relaxation time r. As men- 
tioned above, the difference between r and To comes from 
the momentum dependence of |ak) states. At any rate, 
the full correction da depends only on r due to the addi- 
tional contribution from <5cr 2 (see Fig. [2]: and Eq. (|Bip ). 
Equation (|B1|) can be evaluated using the same proce- 
dure as for Sai . For instance, in Appendix [B] we derive 

{ if (e F -M)/M « 1 , . 
d<T2 ~ \ -(l/3)^i if (e F - M)/M » 1. (A24) 

The full form of the quantum correction, 8a\ + <5cr 2 , de- 
pends only on the transport mean free path r and has (in 
appropriate limits) a universal magnitude, see Eqs. (|2.33[) 
and (j2~35]) . 



C. Calculations 

The road map to 8a starts from a calculation of U in 
Eq. (|2.19p . In Appendix ICl we derive 

LCV = a 8 m i8 mn > + ^ 6 M ^'i' S mi 

+ £ KlSm'V + £ V KnlKn'V , (2-25) 

where n,v G {1,2,3,4}, A ?; = a l r x for i G {1,2,3} 
and A 4 = l 2 r 2 . In addition, a, 6^, c„ and c^,, are 
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Q-dcpendcnt coefficients whose explicit expressions are 
shown in Appendix[C] With those, U is fully determined. 
The next task is to get C™,^,(Q) from Eq. (f2~2Tj) . 

While (lie — U{Q)) can be inverted numerically, it is 
not possible to do so analytically for Q ^ 0. Since we 
are interested in analytical expressions, we follow an ap- 
proximate three-step route. 

First, we diagonalizc (li6 — U) for Q = 0, analytically. 
All eigenvalues can be written in the form A s to, where 
A 9 is the "intrinsic" Cooperon gap or mass. We find that 
one of the eigenvalues has A g = for any e F and M, 
which is a reflection of combined time-reversal symmetry 
and charge conservation. As we elaborate in the next sub- 
section, there may be additional eigenvalues with A g ~ 
when (e F - M)/M < 1 and (e F - M)/M > 1. Here- 
after we refer to eigenvectors of A g ~ eigenvalues as 
gapless (or massless, or "soft") modes. Because A g ~ 
eigenvalues make C large, 5a is determined mainly by 
soft modes. 

Second, we extrapolate the Q = case to Q ^ per- 
turbatively, with the objective of finding how the eigen- 
values of the gapless modes depend on Q. To that end 
6U(Q) = U(0) - U(Q) is written in the basis that di- 
agonalizes U(0). The shift of Q = eigenvalues under 
6U(Q) is then evaluated through standard second order 
perturbation theory. The need to go to second order in 
5U originates from the fact that several matrix elements 
of U™^,(Q) are linear in Q (see Appendix [Cjl . When 
(e F - M)/M < 1 and (e F - M)/M > 1, perturbation 
theory leads to eigenvalues (DQ 2 + A g )rQ. The fact that 
D contains the transport time r rather than the scatter- 
ing time To is generally crucial in order to arrive at the 
correct result for 5a. 

Third, we invert the diagonalized matrix, and trans- 
form its outcome to the \m, m!) basis by using the Q = 
eigenvector matrix (the change of unperturbed eigenvec- 
tors under 6U(Q) is deemed unimportant.) This yields 

Once we have C™,",(Q), we use Eq. (|2.17|) in order to 
extract Cf a (k, k', Q). This is then plugged in Eqs. (j2T22j) 
and (ED. 



D. Results 



term acts somewhat like a Zeeman field in a 2D electron 
gas with Rashba spin-orbit interaction. 

Importantly, the diagonalization of Eq. (|2.21[) reveals 
two qualitatively distinct regimes of quantum interfer- 
ence, (e F - M)/M <C 1 and (e F - M)/M > 1, which 
potentially host additional gapless Cooperon modes. As 
we discuss below, these additional gapless modes can 
change and even reverse the contribution to 5a coming 
from Eq. (gig) . 

When (e F — M)/M 3> 1, we identify a slightly gapped 
(soft) Cooperon mode with 



A g = 2(M 2 /e 2 F )T 



1 = T„ 1 < 7h" 



(2.27) 



whose eigenvector is a spin-singlet and an orbital-triplet: 

±QTB} + \BT))(\U)- lit})- (2.28) 

Physically, r^ 1 is the rate of "intervalley" transitions 
(|T) + \B) -> |T) - \B)) induced by the " mass t erm" 
(Mt z ) in Eq. (pTTj) . Because both Eq. ([2^6]) and 
Eq. (|2.28p are spin-singlets, their contributions to 6a are 
of WAL type (this is proven below). 

Incidentally, M = is the physically relevant regime 
for Weyl scmimctals, which have two degenerate Dirac 
points with linear energy dispersion along the three mo- 
menta axes. Unlike in graphene,— where there are 4 
gapless Cooperon modes (in absence of atomically sharp 
defects and hexagonal warping) , in a Weyl semimetal we 
obtain only 2 gapless Cooperon modes. This difference 
stems from the fact that the SU(2) "valley symmetry" 
of graphene^ gets reduced to a U(l) symmetry in Weyl 
scmimctals, due to the band dispersion along z. Act- 
ing somewhat like a Zeeman field would in a free elec- 
tron gas, the k z dispersion generates a mass for orbital- 
singlet modes, which is why the nearly-gapless Cooperons 
in Eq. (|2~26)) and (f2T28|) are orbital-triplets. 

When (e F — M)/M <t; 1, there are three soft modes 
with gap 



A ff = (2/9)(l-M/e F ) 



2^-1 — 



(2.29) 



Physically, t~ 1 is the rate of spin- flip transitions induced 
by the "spin-orbit term" (i>k • ar x ) in Eq. (|2.1| . The 
eigenvectors for the three slightly gapped modes are 



The diagonalization of Eq. (|2.21[) at Q = shows one 
genuinely gapless Cooperon mode (A g = 0, c.f. Sec. 
IIC), with a spin-singlet and orbital-triplet eigenvector: 



e F + M 



2Ve 



M- 



:\TT) 



e F - M 



M 2 



\BB) 



(in>-ut». 



(2.26) 

The fact that Eq. (|2.26[) remains gapless for any e F /M is 
a physical manifestation of charge conservation. This sit- 
uation differs qualitatively from 2D TIs in HgTc quantum 
wells^i where a nonzero mass term gaps all Cooperons. 
The reason for the difference is that in 2D TIs the mass 



(Ai|2T> + A a |BB))|^) 

(Ai|2T) + A a |BB))|tt> 
(\ 3 \TT) + \ 4 \BB))(\U)- 



lt». 



(2.30) 



where Ai,...,A4 are coefficients that depend only on 
e F /M, such that A x ~ A 3 ~ 1 + 0[(e F /M - l) 2 ] and 
A 2 ~ A 4 ~ 0[(e F /M - 1)]. Therefore, the three soft 
modes in Eq. (|2.30[) are all spin and orbital triplets. As 
will be demonstrated momentarily, their contribution to 
6a is of WL type. 

Next we determine C (c.f. Eq. (|2.23[1 ) by diagonalizing 
Eq. (|2.18[) at Q ^ and doing the angular integration in 
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Eq. (|2~23|) . For (e F - M)/M < 1 we obtain 



67WT 2 



£Q 2 



For (e^ — M)/M » 1, we instead get 



(2.31) 



O 



■■ill 



1 



1 



DQ 2 + t^ l DQ 2 + + Ty 

(2.32) 

In the derivation of Eqs. (|2.31[) and ()2.32j) we have 
included the phase relaxation time and exploited 
DQ 2 t < 1. 

The first term in the square brackets of Eqs. (|2.31[) and 
(|2.32|) is large at Q — > irrespective of € F /M, and origi- 
nates from the spin-singlct Coopcron mode in Eq. (|2.26p . 
Its negative sign means that it makes a contribution to- 
wards WAL. 

Equation (|2.31|) displays a competition between WL 
and WAL, which is no different from that found in an 
ordinary metal with spin-orbit interactions. WL terms 
originate from the three spin triplet modes of Eq. (|2.30[) . 
WL prevails if t7 S> t" 1 , whereas WAL rules if t^ 1 <C 



Equation (|2.32p unveils two different regimes of WAL. 
On one hand, if t7 3> r~ x , the spin-singlet Cooperon 
mode of Eq. (|2.28[) makes a contribution to 8a that equals 
that of Eq. (|2.26|) . In this limit, quantum interference 
can be interpreted as coming from two identical and 
nearly-decoupled Dirac valleys. On the other hand, if 
T7 1 <C Ty 1 , the contribution from Eq. (|2.28l) becomes 
relatively unimportant and the magnitude of WAL is 
halved. In other words, when the intervalley transition 
rate induced by the mass term Mt z is fast compared to 
the phase relaxation rate, the two valleys contribute as 
one. This is quite different from graphene, where strong 
intervalley scattering changes WAL into WL^S The un- 
derlying reason for such a qualitative difference is that in 
graphene a gapless valley-singlet mode is responsible for 
producing WL, whereas in a Weyl semimetal the valley- 
singlct Cooperons are strongly gapped by the k z band 
dispersion. 

Substituting Eqs. ([2~3"Tj) and ([2~5"2"j) in Eq. (|2~2"3|) and 

doing the Q-integral, we arrive at 



6G~aG q lnfo/r) 

-1 if T0 < T s 

a = { 1/2 if r > (t v ,t s ) 

1 if T0<T u , 



(2.33) 



where 5G = WSa is the quantum interference correction 
to conductance and 



G q = e 2 /{2ir 2 h) 



(2.34) 



is a universal conductance unit. In the derivation of 
Eq. ([2331) we have used Eq. ([2^4]) . The reason why 



a = 1/2 when ^> (t v ,t s ) is that in such regime there 
is only one gapless Cooperon mode (hence |a| = 1/2), 
which is a spin-singlet (hence a = \a\). 

While Eq. (|2.33[) is valid in absence of external mag- 
netic fields, the magnetoconductance AG(H) = G(H) — 
G(Q) ~ 5G(H) — 5G(0) can be easily obtained from 
Eq. (j2"lS3"j) for H perpendicular to the TI thin film. The 
replacement of J d 2 Q by an appropriate sum over Landau 
levels^ results in 



AG ~ a GqfiH^/H) 

-1 if Tjj <C T s 
1/2 if T H > (t v ,t s ) 
1 if T H < T v , 



(2.35) 



where f(z) = lnz — "0(1/2 + z) with asymptotes f(z) oc 
z~ 2 for z ^ 1 and f(z) oc ln(l/z) for z -C 1, ip is the 
digamma function, 

t h x = t- 1 + 2eDH/h and = h/ (AcDt^). (2.36) 

Three conclusions of experimental relevance can be ex- 
tracted from Eqs. (|2.33p and (|2.35p . which apply when 
highest occupied electronic states are all located near the 
r point. First, bulk TI bands can display a = — 1 (WL) 
as long as the chemical potential is sufficiently close to 
the bottom of the bulk conduction band. Second, bulk 
TI bands can produce a — 1 when ep/M is sufficiently 
large. Third, when (tp—M^/M is neither large nor small, 
a = 1/2 ensues; this is the conventional result expected 
for ordinary conducting thin films with strong spin-orbit 
coupling, and is the one that has been often presumed 
in experiments on TI films^Ji At th — t s there is a 
crossover between a = — 1 and a = 1/2; likewise, at 
th — t v there is a crossover between a — 1/2 and a = 1. 

The particular expressions for t s and t v in Eqs. (|2.27p 
and (|2.29p rely on our assumption of Vdis oc I4. Spin- 
orbit coupled impurities and/or atomically sharp disor- 
der potentials would induce additional spin- and valley- 
flip processes, whose rates t~j and rj^ 1 would have to be 
incorporated via r^T 1 — > r^T 1 + t~} and r^ 1 -> r" 1 +T~y L . 
If T^ 1 and t~j are strong enough and insensitive to the 
value of ep/M, then the only surviving regime of inter- 
ference corrections is the conventional a = 1/2. 

The conventional a — 1/2 can be found in a wide range 
of parameter space at low temperatures, whereas the un- 
conventional a = — 1 and a = 1 emerge in the rela- 
tively narrow regimes r <C Th t s and r <C th <C t v 
(respectively). How accessible are these unconventional 
regimes? Suppose M ~ 150meV, v ~ 5 x 10 5 m/s and a 
bulk carrier density of n ~ 3 x 10 18 cm~ 3 . This situation 
corresponds to having a small bulk Fermi surface. Then, 
it follows that a ~ — 1 for a fairly wide range of mag- 
netic fields (l H /(12l) < 1, where l H = {Dth) 1 / 2 ). The 
limit a — > 1 is not accesible in this regime. Instead, a ~ 1 
should be accessible in (i) Weyl semimctals or in TIs with 
very narrow bandgaps, (ii) in TIs with large bandgap but 
M(kp) ~ 0. For the latter case it must be kept in mind 



that in the absence of spherical symmetry M(kp) is not 
constant on the Fermi surface. Suppose M ~ 5meV and 
a bulk carrier density of ~ 2 x 10 18 cm -3 . Then, a ~ 1 
in the range of fields corresponding to Ih/{WI) <C 1. For 
typical thin films, the requirements for a = ±1 are com- 
patible with UfI 3> 1. 

Materials like BiTl(Si_iSc5)2, where controlled 
changes of 5 can tune M from to large values,— ap- 
pear to be good candidates to observe crossovers between 
different regimes of magnetoresistance in Eq. (|2.35l) . 

Our analysis has thus far neglected surface states of 
the TI, which can also contribute to the measured mag- 
netoresistance. It can be argued that surface states 
are unimportant and Eq. (|2.35[) suffices in trivial insu- 
lators described by Eq. (|2.1[) . as well as in time-reversal- 
invariant Weyl semimetals and in TIs with very small 
bulk bandgaps (< H/tq). In contrast, when the surface 
states of the TI are robust, Eq. (|2.35[) is incomplete and 
must be generalized. Such generalization is the subject 
for the rest of this paper. 
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FIG. 3: Diagrammatic representation for Soij, where i and j 
are layer indices. For 2D layers without spin-orbit coupling, 
the Cooperon matrix elements are fully characterized by layer 
indices. The velocity operator is diagonal in the layer index; 
therefore, the Cooperons C^k and Off do not enter in the 
expression for Soij. 




FIG. 4: Typical microscopic process that gives rise to 5ai2- 
It can be neglected when the intralayer disorder potentials in 
the two layers are uncorrelated. 



III. QUANTUM CORRECTIONS TO 
CONDUCTIVITY FROM COUPLED BULK AND 
SURFACE STATES 

In this section we consider the combined bulk-surface 
contribution to 5a in 3D TIs with relatively large 
bandgaps. We concentrate on a particular setup that 
consists of a TI thin film gated on one surface. The 
gate voltage can produce a depletion layer that spatially 
separates bulk and surface carriers (Fig. and carriers 
tunnel back and forth across the depletion layer. We as- 
sume the bulk-surface tunneling rate to be much smaller 
than the elastic scattering time on either side of the de- 
pletion layer, so that electrons scatter many times within 
the bulk (surface) before tunneling to the surface (bulk). 
This assumption is experimentally realistic, and it sim- 
plifies the microscopic theory of this section considerably. 



A. Single isolated TI surface 

As a preliminary step, we recall the expression for 5a 
on a single TI surface that is decoupled from the bulk. 
Taking €f s t 3> 1, where cf s is the Fermi energy mea- 
sured from the Dirac point of the surface states, one 
arrive a 16 ' 17 at 

AG/G q = (l/2)f(H <p /H) (3.1) 

for any th- The prefactor 1/2 is consistent with having 
a gapless spin-singlet Cooperon (the spin-triplet Cooper- 
ons have large gaps due to the strong spin-momentum 
coupling on the surface). 



B. Two coupled 2D layers without spin-orbit 
coupling 

As another preliminary step, here we compute 5a for 
two ordinary metallic 2D layers separated by a tunnel 
barrier. In a double layer system, the current flowing 
in layer i can be written as j; = V] ■ a, , E , . where Ej 
is the electric field in layer j. For concreteness we take 
Ei = E2 = E, so that the measured current is j = ji + 
j-2 = cE with a = Y]^ <J%j- Consequently, the quantum 
corrections to conductivity are 5a = . £<xy. The goal 
of this section is to compute 5a from microscopic theory. 

The interference correction 5a^ has the diagrammatic 
representation shown in Fig. [3] Because the velocity op- 
erator is diagonal in the layer index, the only Cooper- 
ons that enter in the conductivity are C^\, with i,j <E 
{1,2}. In particular San involves intralayer Cooperons 
C\\, whereas 5cr 12 and 5a<i\ involve interlayer Cooperons 
C21 and (Fig- HI- Assuming that disorder potentials 
in the two layers are uncorrelated, = for i =/= j. 
This is a reasonable assumption when electrons in the 
two layers scatter off different sets of impurities. Hence, 
we are left with 5a = 5a a. From here on we simplify 
the notation via C\\ = Q. 

When evaluating 5au we will neglect spin-orbit inter- 
actions; however, the main lessons learned in this subsec- 
tion will be transferrable to the spin-orbit coupled case 
studied in the next subsection. In absence of interlayer 
coupling, a standard calculation yields 

^=,-4^A4| Q Ci 0) (Q), (3.2) 

where J Q = J d 2 Q / '(2ir) 2 , an extra factor of 2 is due to 
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(b) 



(a) 



C, 



(0) 



i . i 

i ! i 



C 



(0) 



> 



(b) 



FIG. 5: Single-particle Green's functions, (a) Dressing of 
Bloch states due to intralayer impurity scattering, (b) Dress- 
ing of disorder-averaged Green's functions due to interlayer 
tunneling. The tunneling amplitude is regarded as a random 
variable. 



(0) 



(c) 



i : i 



(0) 



spin degeneracy, Tdi is the scattering time in layer i due 
to elastic impurities (we assume purely s-wave scattering, 
so that there is no difference between the transport scat- 
tering time and the quantum lifetime), is the density 
of states per unit area in layer i and 



cf>(Q) 



1 



(3.3) 



ej>i 



is the Cooperon for an isolated layer. In presence of in- 
terlayer tunneling, C^ in Eq. (|3.2p is replaced by Cf. 



8<Ti, 



Ci(Q), 



(3.4) 



Q 



in whose prefactor we have neglected terms containing 
the ratio between the tunneling rate and the elastic scat- 
tering rate. 

In order to compute Ci, we recognize that the influ- 
ence of interlayer coupling occurs at two different lev- 
els. On one hand, it modifies the single-particle Green's 
function for each layer (Fig. [5]). Because the thickness 
of the depiction layer typically shows microscopic varia- 
tions within the same film as well as from sample to sam- 
ple, the interlayer tunneling amplitude can be regarded 
as a random variable. Consequently, the change in the 
ensemble- averaged Green's function due to tunneling can 



be captured via r 



dl 



-1 



where 



(2n/h)(\t\ 2 )Su j 



(3.5) 



is the tunneling rate from layer i onto layer j ^ i, {\t\ ) 
is the averaged square of the tunneling matrix element 
and 5* is the layer area. Note that (\t\ 2 ) scales like 
so that T t 7 is independent of the layer area. 

On the other hand, interlayer tunneling modifies 
particle-particle correlations that build up Cooperons. 
An approximate diagrammatic expression for these cor- 
relations is shown in Fig. [6] The equation of Fig. [Bfc can 
be solved in momentum space and it yields 



C, 



DjQ 2 



— i 

T 4>3 



{D x Qi+T7l){D 2 Q 



r" 1 ) 



'a 



V 1 



42 

(3.6) 





















c, 






+ 




• 


r(0) 

'"2 


• 


,=r(0) 





















FIG. 6: (a) Cooperon without interlayer tunneling, (b) 
Partially dressed Cooperon Cf^ , where tunneling is included 
solely in the single-particle Green's functions. C 1 } can be 
directly obtained from c} ' via r^i — > f^i. (c) Fully dressed 
Cooperon d, where tunneling is incorporated both in the 
single-particle Green's function and in the particle-particle 
correlations. 




FIG. 7: Typical processes not included in Fig. [6] as they are 
subdominant for tu Tdi- 



for j =/= i. In the derivation of Eq. (|3.6[) we have intro- 
duced 



t- 1 



'<f>i 



(3.7) 



as an effective phase relaxation rate that incorporates 
tunneling, and have used 

/ Gf (k)Gf (-k + Q) c (l-^i- fl,QV dl ) . 

Microscopic processes depicted in Fig. [5] leave out those 
in which two consecutive tunneling events occur without 
any intralayer scattering in between. Likewise, they ig- 
nore electron trajectories in which a tunneling event pre- 
cedes any intralayer scattering (Fig. [7]) . These processes 
are relatively unimportant if r t i 3> Tdi- Not surprisingly, 
Eq. (|3.6[) arises in the coupled equations for the classical 
diffusive conductivity as well (see Appendix iDj) . 
It is convenient to rewrite Ci in Eq. (|3.6p as 



Ci 



B, 



<i 2 



Q 2 



(3. 
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where 




;2 72 
L tl i t2 



(3.9) 



and 
-4, 



1 - B, = {l^ - q 2 a )/(ql - q 2 a ) for j ? i. (3.10) 



In Eq. (|3.9[) we have defined Z^, = (DiT^) 1 / 2 as an ef- 
fective coherence length and lu = (DiTu) 1 ^ 2 as the in- 
tcrlayer leakage length. Besides, q 2 {ql) gets the positive 
(negative) sign in front of the square root. Combining 
Eq. (pT4|) with Eq. ([SIS]) and using A 1 +A 2 =B 1 +B 2 = l, 
we get 



Ee 2 /" 
<5(iii = -2— / 
^ Jq 



1 



1 



?2 



% 2 



(3.11) 



Therefore, the low-field magnetoconductance reads 



Act = Acr « = - G q 



+ / 



1*6 
if 



where 



o(6) = ^9a(fe)/( 4e ) 



(3.12) 



(3.13) 



In the limit of very strong tunneling (th/t^ 0), 
Eq. ([3~T2")l becomes Act ~ -G q f(H b /H), as though there 
was a single layer. In the limit of very weak tunneling 
(rti/r^i — > oo), Act is the sum of contributions from two 
independent films. 

It is helpful to understand the weak and strong cou- 
pling regimes in terms of measurable quantities like the 
interlayer conductance per square, 



gt = {2Tre 2 /H)(\t\ 2 )S^ 2 = a Dl /l 2 t 



(3.14) 



where aoi is the Drude conductivity in layer i. For sim- 
plicity suppose that r^i ~ r^2 = t^. In this case the 
crossover from weak to strong tunneling occurs when 

— -j < 1 (crossover condition), (3.15) 

5t'0 cm &D2 

namely when the tunneling resistance for a square of 
area l 2 ^ becomes smaller than the sum of the classical 
intralayer resistivities. Let us define 



9c 



= ( a Dl + <j d\) 1 1 



(3.16) 



If g t < g c , then Act/G 9 ~ -21n(ff/ffy) for H > H+. 
If g t > 5c, then Aa/G q ~ - \n{H/H^) for ffy < if < 
H^{g t /g c )- Thus changing the interlayer conductance re- 
sults in a factor-of-two change for the magnitude of the 
WL correction. 

Limits reminiscent of the above were first discussed in 
inversion layers of multivallcy semiconductors like Sij 2 ^ 



FIG. 8: Example of an interlayer scattering process that is 
allowed in multilayer systems. Its analog in multivalley semi- 
conductors of Ref. [24j is forbidden. 



where the role of the layers is played by different electron 
pockets in the Brillouin zone. Similarities notwithstand- 
ing, there are clear differences between our microscopic 
theory and that of multivalley semiconductors. On one 
hand, the separation in momentum between valleys of Si 
prevents scattering processes such as the one in Fig. [5] 
These processes arc not only allowed in our case, but 
also lead to the Cooperon dressing shown in Fig |BJ:. On 
the other hand, in our case the interlayer Cooperon van- 
ishes due to uncorrelatcd disorder potentials in the two 
spatially separated layers. That is not the case in multi- 
valley semiconductors, where both valleys scatter off the 
same set of real-space impurities and intervalley Cooper- 
ons contribute crucially to 5a. 

Finally, it should be mentioned that Eqs. ([331) . (|3.10j) 
and ()3.12j) coincide with those derived by G. Bergmann^ 
who invoked macroscopic arguments based on the diffu- 
sion equation. The microscopic theory of this subsection 
supports Bcrgmann's results, insofar as tu 3> Tdi and the 
disorder potentials in the two layers are uncorrelated. In- 
cidentally, yet another way to arrive at the same results 
is unveiled in Appendix [Ej this later method will prove 
convenient in the upcoming subsection. 



C. 3D TI film with bulk-surface coupling 

We now consider a 3D TI film (Fig. [S]) with a gate 
electrode placed near its top surface. 

At the moment we neglect the bottom surface of the 
TI, which will be incorporated below. For ease of nota- 
tion we use subscript "1" to refer to "bulk", and sub- 
script "2" to refer to "top surface". Like in the preced- 
ing subsection we assume bulk-surface disorder correla- 
tions to be negligible, so that the quantum corrections 
to conductance can be written as SG = SGn + SG22 = 
WSan + (5<722 • SG is approximately independent of the 
film thickness W as long as W <C l^i , where l^i was 
defined below Eq. (|3T0j) . 

The goal of this subsection is to calculate SG from 
microscopic theory. Unlike in the previous subsection, 
here both "layers" are spin-orbit coupled. We assume 
that tunneling events, albeit being time-reversal invari- 
ant, conserve neither spin nor orbital degrees of freedom. 
Indeed, in a TI spin is not conserved for non-momentum- 
conserving tunneling. Similarly, the orbital degree of 
freedom is not conserved due to broken inversion sym- 
metry near the surface. 
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Let us begin with no tunneling. On one hand, there are 
four surface Cooperon modes: one gapless spin-singlet 
mode and three spin-triplet modes with large (~ 
gaps. On the other hand, there are sixteen bulk Cooper- 
ons, of which a spin-singlet mode (Eq. (|2.26[1 ) is always 
gapless. In addition, four of the bulk modes (the spin- 
singlet of Eq. (fiT28j) and three spin-triplets of Eq. ([230]) ) 
can be "soft" depending on ep/M. The rest of the bulk 
Cooperon modes have large gaps of order t^ 1 . 

Let us now turn on tunneling. Since Tu 3> r^, we can 
limit ourselves to analyzing the effects of tunneling within 
the low-energy subspace formed by the soft Cooperons. 
If there are no magnetic impurities in the depletion layer, 
the total spin of the Cooperon is a good quantum num- 
ber even in presence of tunneling. Accordingly tunneling 
does not mix spin-singlet modes with spin-triplet modes, 
and the full (dressed) Cooperons can also be classified 
into spin-singlets and a spin-triplet. 

In the regimes t^i <C t s and t^i 3> (t s ,t„), tunneling 
dresses one soft spin-singlet Cooperon in the bulk with 
another soft spin-singlet Cooperon on the surface. This 
dressing is completed as explained in Section IIIB: first 
by renormalizing the phase relaxation time — > r^j, 
and afterwards proceeding with the series expansion of 
Fig. [6}:. All "blocks" appearing in this series expansion 
are spin-singlets. When r^i <C t s , the soft spin-triplet 
Cooperons from the bulk are dressed simply through 
7>i —t f^i: they do not get appreciably admixed with 
the spin-triplet Cooperon on the surface because the lat- 
ter has a large gap. 

In the regime t^i <C t v , there are two gapless singlet 
Cooperons in the bulk, each of which can hybridize with 
the singlet gapless Cooperon on the surface. For this 
situation, Fig. |Bt does not capture all possible processes 
and the calculation from the previous subsection must 
be generalized; this generalization is carried out in Ap- 
pendix |e] 

With the above considerations in mind, we combine 
Eqs. (|2.35j) and (|3.1j) in order to obtain the total contri- 
bution to low-field magnetoconductance: 




/(*)+/($)" 3/ (¥) if 

f{it)+f{%) if f H »(T s ,r 1J ) 

/(f)+/(f L )+/(% i ) if f*<r„, 



(3.17) 



where Hi = hqf/(4e) for I G {a, b, c, d}, 



H<f,i = ft/(4e.Dif0i), and t h 1 



2eD 1 H/h. 

(3.18) 

Note that the effective phase relaxation rate increases 
linearly with the bulk-to-surface tunneling rate (c.f. 
Eq. (|3.7[1 ). The characteristic momenta q a (b) have been 
introduced earlier in Eq. p.9[) . The additional momenta 
q c (d) arc identical to q a n,\, except for t^ 1 —> 2r t 2 . The 
reason for this difference is that the surface Cooperon can 
decay into two gapless bulk Cooperons when r^i <C t v . 



The first line of Eq. (|3.17p displays a competition be- 
tween WL and WAL, and suggests that it is possible to 
induce a WAL-to-WL transition with a varying gate volt- 
age. In the weak tunneling regime WL prevails, whereas 
in the strong tunneling regime WAL takes over. Simi- 
larly, a gate voltage can induce transitions between three 
different WAL coefficients: a G (1/2,1) in the second 
line, and a G (1/2, 3/2) in the third line. The second line 
of Eq. (|3.17j) describes quantum corrections as if they 
originated from two independent thin films with mixed 
bulk-surface character; indeed, universal results expected 
for the simplectic symmetry class are recovered when 
the effective phase relaxation times become the longest 
timescales of the problem. Some simple limiting cases of 
Eq. (|3.17|) are discussed in Appendix |FJ 

Thus far we have considered the coupling between the 
bulk and one (the top) surface of the TI film. As a conse- 
quence, Eq. p,17[) applies to a TI film only if the phase 
relaxation time of the bottom surface (adjacent to the 
substrate) is short compared to other phase relaxation 
and tunneling times in the problem. This condition is 
likely not met in some recent experiments^^ which re- 
port on independent contributions from both surfaces. 
Partly motivated by these experiments, we now general- 
ize Eq. p,17[) so as to capture two surfaces, each coupled 
to bulk states. 

We consider the scenario depicted in Fig. |H1 where 
the bottom surface contains bulk carriers. Since there 
is no depletion layer near z = W, we assume that the 
bulk-surface tunneling rate therein is strong compared 
to the phase relaxation rate, yet weak compared to dis- 
order scattering rate. Hence we describe the hybrid of 
bottom surface and bulk states via Eq. (|F1[) . and there- 
after couple this hybrid to the top surface along the lines 
of Eq. (|3.17p . The resulting expression for AG can be 
approximated as 



G q 2 J \ H ) 2 J \ H 



(3.19) 



where H' a rj,\ = ^(<?l(t,)) 2 /(4e)- The characteristic mo- 
menta q' a and q' b obey Eq. (|3.9p . where "1" labels the 
top surface and "2" labels a hybrid between the bottom 
surface and the bulk. 

Notably, Eq. (|3.19[) implies that WL is no longer pos- 
sible once the bottom surface is strongly coupled to 
bulk states. Instead, conventional WAL ensues with 
a G (1/2,1). This observation not only sheds light on 
why current experiments see no indication for WL, but 
it also gives insight as to how WL could be observed in 
TI films. 

A possible strategy is to degrade the surfaces, e.g. by 
depositing magnetic impurities on them, and decoupling 
them from the bulk by double-sided gating. One may 
expect WL even if only the top surface is decoupled, 
while the (degraded) bottom surface is in contact with 
the bulk. In this case, Eq. (|3~T7) reduces to Eq. ([235)) 
derived for the sole bulk conduction, with the replace- 
ment tT^ — > r7y +TJ3 1 , where Ttz is the tunneling rate of 
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FIG. 9: Schematic energy band profile for a gated 3D TI thin 
film, z = corresponds to the top surface of the device, im- 
mediately under the gate, z — W corresponds to the bottom 
(ungated) surface. The vertical (blue) solid lines at z = 0, W 
are surface states. The curved solid (red) line is the bulk con- 
duction band, and the dot-dashed (brown) curve is the bulk 
valence band. The chemical potential is depicted by a hori- 
zontal dashed line. Zd is the thickness of the depletion layer, 
where neither bulk nor surface carriers are present. tp s is the 
Fermi energy of the surface states measured from the Dirac 
point (£f s < in this figure). e_F is the Fermi energy of the 
bulk states, measured with respect to the midgap point. 



electrons from bulk to the bottom surface. If the film is 
thick enough, then may become sufficiently small to 
provide some dynamic range for observing WL behavior. 
This same strategy can also facilitate the observation of 
WAL with a > 1. 



D. Estimates for the bulk-surface coupling 

This subsection is devoted to an approximate electro- 
static and quantum mechanical analysis of the depletion 
layer in a TI film, which will result in quantitative esti- 
mates for the bulk-surface coupling. 

For a TI with an n-doped bulk, a negative charge per 
unit area (—Q g ) placed at the gate repels electrons from 
bulk bands at z = as well as from the surface states 
at z = 0. This leaves a positive net charge on the top 
surface, which is equivalent to a downward shift in the 
local chemical potential at z = 0: A// s = ep — ep s . Re- 
call that cfs is the Fermi energy of the surface states 
measured from the Dirac point (for simplicity the Dirac 
point is assumed to be in the middle of the bandgap at 
z = 0) and ep is the Fermi energy of the bulk states mea- 
sured from the middle of the bandgap. Since the chemi- 
cal potential deep inside the bulk must be unaffected by 
the gate, Afj, s ^ implies a band bending of magnitude 
(f> s = A/z s near the gated surface (Fig. [SJ. 

When A/i s > (ep — M) there are no bulk carriers left 
at z = and a depletion layer appears at z € (0,Zd), 
where Zd is determined below. For each value of Q gi A/z s 
(or equivalently ep s ) can be uniquely determined from 
the overall neutrality condition Q s + Qd = Q g , where Q s 
is the positive net charge induced on the surface, and Qd 
is the positive net charge in the depletion layer. 



where n (c.f. Eq. I|2.3[) ) is equal to the density of charged 
donors in the depleted region. The electrostatic energy 
profile in the depleted region then obeys 

1 P 71 

<$>{z) = 4> h -- — (z-zd)\ (3.20) 

where = 4> s — {ep — M) = M — ep s , k is the static 
dielectric constant and 



Zd 



(3.21) 



In the derivation of Eq. (|3.20p we have assumed that 
the electric field vanishes at z = Zd, which is accurate 
within a screening radius. As the gate voltage is made 
more negative, the maximum width of the depletion layer 
(z™ ax ) is achieved when <j) b ~ 2M. For (j) b > 2M, the 
bulk bands get inverted at z = and Zd saturates. We 
estimate z™ ax ~ 20nm for some typical parameter values 



(M = 150meV, n ~ 4 x 10 18 cm 



50). 



Once the electrostatic profile of the TI film is charac- 
terized, we can analyze the quantum mechanical tunnel- 
ing of electrons across the depletion layer. The tunneling 
conductance per unit area is roughly 



9i 



(e 2 /h)X- 2 exp(-2 X ), 



(3.22) 



where Xp is the smallest between bulk and surface Fermi 
wavelengths, and 



X 



dz- 



hv 



1 e 2 nz\ 
6 k Hv 



(3.23) 



In Eq. (|3.23[) we have ignored effective mass and Fermi ve- 
locity mismatches across the depletion layer. The WKB 
exponent % can be tuned by a gate voltage: as Zd varies 
from to z™ ax , \ goes from to ~ 6. 

Drawing from the previous subsection (c.f. Eq. (|3.15jl ). 
the crossover from weak to strong bulk-surface coupling 
occurs when 



1 

~9A 



< 



1 



1 



(TDlW (JD2 



1 

0"_D2 



~ , (3.24) 



where in the second equality we have assumed that 
ajjiW 3> &D2- This is a good assumption provided that 
(i) the bulk mean free path is of the same order as the 
surface mean free path, and (ii) kpW ^> 1. Plugging 
Eq. (|3~2"2"|) in Eq. QSTflfy , the latter becomes 



^>(W 2 ) 1/2 cxp( X ), 

Ap 



(3.25) 



In the depletion approximation^* one has Q 



nz d , 



where kp s = \ep s \/hv is the Fermi wave vector for the 
surface states and we have used ao2 ~ (e 2 /h)kp s l2- 

When z d = z£ lax , the right hand side of Eq. (|3~251> 
reaches ~ 1000, which exceeds the typical I^/Xf in TI 
thin films by at least an order of magnitude. Therefore, 
when the depletion layer has its maximum width, the 
top surface and the bulk of the TI film can be regarded 
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as weakly coupled. This state of affairs changes rapidly 
when the depletion layer is made thinner by a gate volt- 
age. For instance, when z c i = z^^/y/ 7 !, the right hand 
side of Eq. (|3.25p equals ~ 30, which is comparable to the 
typical l^/Xp. Further slight reductions in Zd can subse- 
quently drive the film into a regime of strong bulk-surface 
coupling. These estimates justify the interpretation of 
experimental data given in e.g. Ref. fill ]. 



IV. SUMMARY AND CONCLUSIONS 

We have completed a theoretical study of low- 
field magnetoresistance in electrostatically gated 3D TI 
films. The concise analytical expressions presented here 
[Eqs. ([235]) . (|3"T7|) and (j3TT§)) ] may shed light on the 
quantum magnetoresistance of TIs, Weyl semimetals, as 
well as some topologically trivial materials. Only mag- 
netic fields that are perpendicular to the TI thin film 
have been considered in this work; for in-plane fields 
and small bulk bandgaps, quantum interference contri- 
butions might be masked by classical magnetoresistance 
anomalies i^i 

A number of predictions from this work have not been 
articulated in previous studies and await experimental 
confirmation. For instance, we find that TI thin films 
with low bulk doping may exhibit weak localization (WL) 
or negative magnetoresistance, instead of the often pre- 
sumed weak antilocalization (WAL) or positive magne- 
toresistance. Admittedly, the parameter space for WL 
is relatively narrow, and vanishes when either surface of 



J 



the TI film is strongly coupled to bulk states. However, 
WL may be experimentally accessible in thicker films, or 
in thin films where the surfaces have short phase relax- 
ation times. Under these conditions, a gate can induce 
a crossover between WL and WAL. On a separate note, 
we find that the "universal" prefactor for WAL varies de- 
pending on the bandgap of the TI, on the bulk doping 
concentration, on the phase relaxation times, and on the 
applied gate voltage. 

The results from this work are applicable to conduct- 
ing yet lighly doped TIs, with thicknesses ranging be- 
tween the bulk transport mean free path and the bulk 
phase relaxation length. It may be useful to find out 
how the results derived here change in highly doped TIs 
containing additional electrons pockets away from the T 
point. Likewise, it may be helpful to extend our results to 
thinner films. Other potentially interesting tasks involve 
investigating universal conductance fluctuations and de- 
termining the influence of electron-electron interactions 
in the magnetoresistance of doped TI films. 
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Appendix A: Renormalized velocity operator 

The velocity operators appearing in the expressions for ud and 5a (c.f. Sec. IIB) must be renormalized with ladder 
diagrams containing impurity scattering. The Dyson equation for the renormalized velocity operator is (Fig. 1 10[) 



v Q/J (k)=v Q/3 (k)+u ]T / (ak\a'k')(f3'k'\(3k)G A (k')G R (k')v a , p ,(k'), 



(Al) 



where v Q( g(k) = S a phv 2 k/Ek is a matrix element for the bare velocity operator. We solve Eq. (|A1[) by guessing a 
solution of the form 



v Q/3 (k) = 7/ck<5 Q/ j, 



(A2) 



where jh is a scalar that depends on |k| but not k. Although it is a priori not obvious that the renormalized velocity 
operator should be diagonal in the band indices, substituting Eq. (| A2|) in Eq. (|A1|) and using Eq. (|2.4[) we find that 
v Q( g(k) oc 5 a p is indeed appropriate provided that 



Ik = • 

Eh T Q 



Here 



(A3) 



*=2«vu [ ^l^|( akF |a'k' F )| 2 (l-k F .k^) (A4) 

a.' 

is the transport scattering time. Therefore, the final result for the renormalized velocity is v a p(k) = v a/ 3(k)(r/ro). 
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FIG. 10: Impurity vertex corrections for the velocity operator 



Appendix B: Evaluation of <5o2 in some simple cases 
The expression for 5<72 (depicted in Fig. [2J:) reads 



5a 2 c-2^[ ^(k)^(kOG^(k)G^(kOG^(-k)G^(-k')G«(-kOG fl (k) £ r^(k,-k',0)l /0C a %(k,k', 

(Bl) 

where the overall factor of two stems from the fact that the two diagrams in Fig. [2t give identical contribution, and 
the band indices a, (3 etc. are summed over 1, 2. For generic (ep — M) /M, the calculation of 5a 2 is cumbersome. Here 
we focus on two simple limits that are of interest: (ep — M)/M -C 1 and (ep — M)/M 3> 1. 

When (ep — M)/M <C 1, the momentum dependence of |akp) is negligible. Consequently, C^,(kp, k' F , Q) and 
r^fg (kp, — k' F , 0) become independent of k^? and k^.. Since the matrix elements of the velocity operator are odd 
under k — > k and k' — > — k', it is clear that 



<5ct 2 ~ 0. 



(B2) 



The limit of (ep — M) /M 3> 1 is less trivial. In this regime the Hamiltonian is approximately block diagonal both in 
absence and in presence of disorder, because the disorder potential we take is spin- and orbital-indpcndcnt. Therefore 
we may focus on a 2 x 2 Hamiltonian (describing a Weyl node of positive chirality), 



h'(k) =fek-<T + Vo(r)l 2x2 , 



(B3) 



where k = fc(sin 9 cos 4>, sin 9 sine/), cos 9). The result for 8a obtained from such Hamiltonian needs to be multi- 
plied by two at the end, as each block makes an equal contribution. The eigenstates for h'(k) are |+,k) = 
(cos(6»/2),exp(^)sin(6'/2)) T and |-,k) = (sin(0/2), - exp(i0) cos((9/2)) T . 

One significant simplification from Eq. (|B3[) is that there is only one band at the Fermi energy. This allows us to 
rewrite Eq. (|B1[) as 



5 a 2 



2tt 
dVL 



-in 



T 2 




Uq „ 


[/ 






k f dn k , 


r J 4tt 



dkk 2 kv 2 



(G A ) 2 G R 



1 

W 



d 2 Q 
(2n) 2 



sin 9 cos <t> sin 9' cos </>'(+, k F \+, k' F ) (+, -k F \+, -k' F )Cll(k F , k' F , Q), 



(B4) 



where the aforementioned extra factor of two has been accounted for. It is illustrative to compare Eq. (|B4j) with its 
counterpart in 5a\: 



dkk 2 k 2 v 4 



1 



d 2 Q 



(B5) 



In Section II we detailed the steps to follow for the evaluation of Eq. (|B5I) . Applying those same steps to Eq. (|B4 
and using 



dkk 2 k 2 
^E 2 



R\2tr<A\2 



(G ) (G 



4tt^t 3 

Try and 



(ifcfc 2 fc 



(G A ) 2 G R 



we arrive at 



(2tt) 2 £ fe 
<5er 2 = "g^ ! = "3^9 m (~) 



u 



h 5 v 2 



(B6) 



(B7) 
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Appendix C: Evaluation of matrix elements for U 

In this Appendix we calculate the coefficients entering in Eq. (|2 . 25|) . These coefficients generally depend on the 
frequency f2 and wave vector Q of the external perturbation. Even though only Q — is needed for our evaluation of 
5a, for completeness here we allow for f2 ^ as well. 

The calculation is facilitated by rewriting Eq. (|2.1[) as 



ft(k) = J>„(k)A' 



(CI) 



where /it g {1,2,3,4}, ^(k) = hvk, and A 1 = a l r x for i g {1,2,3}, 7j 4 (k) = M and A 4 = 1 2 t z . Then, the 
finite-frequency retarded and advanced Green's functions read 



G£ ( „ A) (k,fi) 



[e*M]*-El 



(C2) 



where e fl = e F + 17 and e A ee e F + Ml - 17, with 7 = ft/(2r ) (c.f. Eq. ([2~5j) 7 Substituting Eq. ([U2j in Eq. (|2~T9|) . we 

get 



(C3) 



where 



a = ii 



c M = mo 



d 3 k 



e R (e A + Ml) 



f d 3 k 

; b » = u °J 



e%(k + Q) 



)3 [ {e xy-E\][(eA + hny-E2 +Q ] 



d 3 k 



(e A + Ml)d fl (-k) 



(27T)3 [( e «)2_ jB 2 k][(e A + m) 2_ jB 2 +Q] 



li 



d M (-k)d v (k + Q) 



(27T)3 [ (e «)2_ £ ;2 k][( ^ + ^ )2 _ jB 2 +Qr 



(C4) 



and /i, ^ G {1, 2, 3, 4}. In the diffusive transport regime, namely (ep — M) ^> 7 ^> (hvQ, Ml), the integrals in Eq. (|C4 
can be analytically performed and the outcome is 



,.(0) 



1-1(1 
12 



ai 2 \ h 2 v 2 Q 2 inn 

~ 27 



h = -ci ~ 1 



£ (°) - —\ ^ v Q x 



dxi 



d 22 * ~ I 
-ill 



A/ 5 



M" 



M 



«33 

.5 

A/ 2 

"44 — — t - a 

ei 



,(°) 



,(o) 



,(°) 



7 
1 

1 

20 

1-^-11 
20 



-c 2 ~ 


J.<«(l- 


A/ 2 \ 


(3Q 2 


+ Ql)h 2 v 2 


iMl 




7 2 






f Q 2 )ftV 


iMl 




7 2 





M 2 \ /1VQ 2 iMl 



^12 = ^21 — a 



(0). 



30 



M 2 \ h 2 v 2 Q x Q 



di4 = —da ~ — ia^ — ( 1 5- 



r 

A/ 2 \ ftwQ^ 



7 M 

04 = C4 = — a 



7 



-^42 ~ — ior 



(Q) M ( 1 M 2 \hvQ y 



7 



(C5) 



where = [2(1 + M 2 /ef.)] _1 , and the elements omitted above are zero. It is worth noting that Eq. (|C5[) can be 
used to investigate the dynamical spin-charge coupling in doped TIs. Since this task is not directly related to the 
theme of this paper, it will be pursued elsewhere. 
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Appendix D: Classical conductivity of two coupled layers 

In this Appendix wc analyze the classical conductivity of two coupled layers. The current in layer i is given by 
ji = ~^2j o'ij'Eij- It is illustrative to write er^ in terms of the diffusive density-density response, using the continuity 
equation 

|+V-j + A^(ft-ft)=0 (Dl) 

3 

along with the constitutive equation j; = — D^Vpi — e 2 ViDi~Ei. A is the interlayer tunneling rate. Thus it follows that 

aij(q, u) = ~%Xij + ^2^2(Xij - Xkj), (D2) 

q q k 

where Xy(q,k>) = e VjDjq pij(q,u) is the density-density response function and 

^ (q ' w) -\A^ 0) /(i-A^ o) ) if i*r m 

with pf^ = (D{q 2 — iuj + A) -1 . The dressed diffusion probability pa, derived here from the continuity equation, has 
identical form as Eq. Q3.6p . which was derived microscopically in Section IIIB. Here oj and q are the frequency and 
momentum associated with the applied electric field. A straightforward calculation shows that a\ 2 = a 21 = when 
Ej is spatially uniform (q = 0) . 



Appendix E: Equations for coupled Cooperons 

In the first part of this Appendix we present an alternative derivation for the results of Section IIIB. In the second 
part of the Appendix we generalize the derivation to make it suitable for TI thin films with r^i <C t v , which contain 
two gaplcss singlet Cooperons in the bulk and one gapless singlet Cooperon on the surface. The outcome of such 
generalization is the third line of Eq. (|3.17p . 

1. Two 2D layers without spin-orbit coupling 

In this subsection we use "1" and "2" to label the two layers. The relevant Cooperon modes are then Cu, C12, 
C21 and 622- Recognizing that Cooperons must obey a diffusion equation in absence of phase relaxation, we posit 
the following coupled equations: 

{D1Q 2 + T^)C n + MCn - C21) = V^wJx) 

(D2Q 2 + t^)C 2 i + X(C 21 - C11) = 

(D 2 Q 2 + r^)Qa + MC22 - C 12 ) = h/(2^2T 2 d2 ) 

(D1Q 2 + t^)C 12 + X(C 12 - C22) = 0, (El) 

where A is the interlayer tunneling rate. Note that the source term appears only for the diagonal terms of the 2x2 
Cooperon matrix. The solution of Eq. (|E1[) reads 



C'22 



D 2 Q 2 + f^ 



2nu 1 r 2 {D1Q 2 + f^)(D 2 Q 2 + f^) - X 2 

__h DiQ 2 + r-l 

2^ 2 T 2 2 (D^ 2 + r^)(D 2 Q 2 + f£) - X 2 

D 2 Q 2 +fr' 



C12 — C21 — — ^ ,_ 1 Cu, (E2) 
T 4>2 
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where t^ 1 = r^ 1 + A. The expressions for C\\ and C 22 agree with Eq. (|3.6[) . In addition, C\% and C21 agree with 
the expressions for p\2 and p 2 i derived in Appendix [D] (where we discussed the classical diffusive conductivity). Cu 
of Eq. ((EI) is equivalent to C| of Fig. [3] Likewise, Ci 2 and C 2 i of Eq. (|ET]) correspond to C\\ and C 22 of Fig. ED 
Although C12 and C21 are nonzero, they do not contribute to 8a because the velocity operator is diagonal in the layer 
index. Therefore, we reproduce the expression of Section IIIB for 8a. 



2. TI film with two gapless bulk Cooperons and one gapless surface Cooperon 



In this subsection we use "1" and "3" to label the two bulk Cooperons, and "2" to label the surface Cooperon. The 
generalization of Eq. (|E1[) is 



(D 1 Q 2 +t^)C 11 + X(C 11 -C 2 i) 
(D2Q 2 + T^)Cn + A(2C 21 - C u - Cat) 
{D 1 Q 2 +t^)C 31 + X{C 31 -C 2X ) 



h/{2-Kv ir 2 dl ) 



0. 



(E3) 



(DiQ 2 + r~l)C X 2 + A(C 12 - C22) = 



(D2Q 2 + t^)C 2 2 + A(2C 22 - C 12 - C 



2,2 1 



{D 1 Q 2 +r^)C 32 + X(C 32 -Q 



22 1 



h/{2TTV 2 T 2 2 ) 





(E4) 



and 



(AQ 2 + t^)C 13 + X(C 13 - C 23 ) = 
(D 2 Q 2 + t^)C 23 + X{2C 23 - C13 - C 33 ) = 

(p x Q 2 + r^)C 33 + X(C 33 - C 23 ) = h/(27r^ 1 T 2 1 ), (E5) 

Once again in Eqs. (|E3|) - (|E5|) the source term appears for the diagonal components of the 3x3 Cooperon matrix. 
In addition, a factor of 2 has been multiplied in front of some tunneling rates associated to surface Cooperons. The 
rationale behind this is that the Cooperon on the surface can decay into two bulk modes, i.e. the effective decay rate 



becomes 



rj>2 



2X. Aside from this, we have assumed a unique tunneling rate A between all pairs of Cooperons. 



The quantum correction to conductance can be written as 

/Q "■- JQ 

Solving Eqs. (|E3j) - ()E5j) requires some algebra. The results for the Cooperons of interest are 

h 



8G = 2^u 1 D 1 t 2 1 J (Cn+C 33 )+2^is2D2T 2 2 J C 22 . 



(E6) 



C11 — C33 — 



C22 — 



(D1Q 2 + f^)(D 2 Q 2 + t£ + A) - A 2 



2™!^! {D X Q 2 + f - 1 ) (D.Q 2 + f^)(D 2 Q 2 + f7 2 x + A) - 2A 2 



DiQ 2 



2™ 2 r 2 2 (D1Q 2 + f£)(D 3 Q* + f£ + A) - 2A 2 ' 
which are not illuminating expressions. It is better to rewrite them as 



(E7) 



Cu — C; 



1 



:!.'! 



X 



Y 



Z 



r + £ 



r+q 2 



C22 

so that Eq. (|E6I) transforms into 

e 2 

SG = — 



1 



27w 2 t 2 2 D 2 



A 



B 



Q 



X 



Q 2 +ql Q 2 



Y 



z 



A 



D 



f + q 2 



r + q 2 



<P + q* Q 2 + q 2 a Q 



(E8) 



(E9) 
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Comparing Eqs. (|E7[) and (|E8|). we arrive at 

_ i 



ql 



1 9 



X = 



Y 



(Diq 2 x 



V 1 



; s = i 



A 



A) -A 2 



A#2(g 2 -<z 2 )(<7 2 -g 2 ) 

=-1/^-1 , \\ , n „ 



2--1 



Vl 1 ^ + A ) + Diq 2 y (-D 2 q 2 y + t£ + A) + A 2 



i?!^ 2 -^ 2 -^ 2 ) 

^x(rJ + A) + ^ig 2 (-^2^ 2 + fJa 1 + A) + A 



D 1 D 2 {ql-ql){ql-ql) 



(E10) 



and 



A 



(6) 



2 _ 2 



q 2 z = l/(£>!^ 



± 



A 

£2" 



8A 2 

^1^2 



(Ell) 

Note that q a (b) = qx(y)i which will be important below. Also note that the expressions for A, B and q a (b) are identical 
to the ones in Section IIIB, except for the following difference: the effective inelastic scattering rate for layer 2 is now 
T 4>2 + ^A instead of t^ 2 + A, for the reason explained above. 

Although Eqs. (|E10[) and (|Elip look cumbersome, after substituting Eq. (|E11[) back in Eq. (|E10[) wc find some 
remarkable simplifications. In particular 



Z = l/2 , 2X + A = 1 and 2Y 
Replacing these in Eq. (|E9[) immediately leads to 



B = 1. 



6G = 



1 



1 



«2 



Q 2 + <z 2 Q 2 + <z? 



In consequence, we recover the third line of Eq. (|3.17p for the low-field magnetoconductance: 



AG 

G„ 



f 



Ha 

H 



+ .f 



Hb 
H 



+ f 



~~H 



where H a = hq 2 /(4e), etc. As a reality check, let us take some simple limits. 

First, consider the case of no bulk-surface coupling, A — > 0. In this case H a 
fi/(4e£'2T'0 2 ), which produces 



AG 

G„ 



H 



2/ -FT + / "FT 



H, 



H 



(E12) 
(E13) 

(E14) 

S/(4e2?iT0i) and Hb = 
(E15) 



This is indeed the result one would have expected when bulk and surface are decoupled. 

Second, suppose both t^i and r^2 are infinitey large, for arbitrary tunneling rate. Then it follows that Hb = 0, 



4e V D x D 2 



Then, 



AG 

G a 



H„ 



H 



and H z = — 

4eDi 



H, 



f "FT +/ 



H 



(E16) 



(E17) 



The fact that Hb = means that wc recover the conventional WAL case (as we should when the phase relaxation 
times are infinitely long). 

Finally, consider the case of very strong tunneling between bulk and surface states. In this case H a and H z become 
very large (oc A), whereas Hb becomes independent of A. Consequently 



AG = 1 /Hb 
2 1 \ H 



G 



(E18) 



as if we had a single channel contributing to WAL. This seems to make sense too, because when tunneling is strong, 
Ca are strongly coupled to one another (i = 1, 2, 3). 
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Appendix F: Some special cases of Eq. (|3.17[) 

In this Appendix we analyze some simple limiting cases of Eq. (|3 . 1 T[) , which considers a single TI surface coupled 
to bulk states. First, suppose that surface-bulk tunneling is strong, so that T t i <C tm for i = 1,2. In this case 
(H a ,H c ,Hi) > (H b ,H d ) and thus Eq. (pTT?|) turns into 

A r 1 f f( H b/ H ) if f H < T S 

=^ = -l f(H b /H) if f H » (t v ,t s ) (Fl) 
Uq 1 { f(H d /H) if th < r v , 

where H b ~ ft/(4e)(l/r i + l/r 02 )/(£»i + D 2 ) and if d ~ ft/(4e)(2/r i + l/T <t>2 )/{2D l + D 2 ). For simplicity we have 
taken T t \ — r t 2, but this assumption can be easily relaxed. In sum, WL is not possible when the bulk-surface coupling 
is strong, and the film exhibits conventional WAL (a — 1/2) regardless of the bulk carrier concentration. 

Next, we consider a weak surface-bulk tunneling, so that tu 3> t^ for i = 1,2. In this case the outcome depends on 
whether D\t^\ > Dit^ 2 or D\t$\ < D%t^. Without loss of generality suppose that D\t$\ > D 2 t^,2- Then Eq. (|3.17[) 
yields 

Ar i { f(H^ 2 /H) - 2f(Ha/H) if th « r s 

— ~ - I m^/H) + f(H^/H) if th » (t v ,t s ) (F2) 
Uq 1 { fiH^/H) + 2f{H^JH) if t h <t v , 

where iJ^i = h/ (AeDiT^i) for i = 1,2. When if^i and i?^2 are of the same order, the first line of Eq. (|F2[) displays 
WL with a = — 1/2 and the third line exhibits WAL with a = 3/2. If instead <C #02, AG is the same as if there 
were no surface states. This latter regime can be experimentally accessible by e.g. depositing magnetic impurities on 
the surface of the TI. 

Last, we consider the case T t \ 3> t^i 3> r t 2 for i = 1,2. This situation may be relevant for some thicker TI films 
where r t i/r t 2 = Wvi/^ ^> 1 (for thicker films, surface states have more bulk states to decay onto). The resulting 
magnctoconductance is once again as though there were no surface states: 

Ac { -f(H^/H) if f H < T a 

— = { \f(H^/H) if f H » (r„,7i) (F3) 
ljq [ f(H<f>i/H) if t h <^t v . 
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